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ABSTRACT 


A model for simulating the measured radar backscatter ing 
coefficient of vegetation-covered soil surfaces is presented 
in this study. The model includes both coherent and inco- 
herent components of the backscattered radar pulses from a 
rough soil surface. This surface is characterized by two 
parameters, the surface height standard deviation a, and 
the horizontal correlation length The effect of vege- 
tation canopy scattering is also incorporated into the model 
by making the radar pulse subject to two-way attenuation and 
volume scattering when it passes through the vegetation 
layer. These processes are characterized by the two param- 
eters, the canopy optical thickness t and the volume 
scattering factor n. The model results agree well with 
the measured angular distributions of the radar backscatter- 
ing coefficient for HH polarization at the 1.6 GHz and 
4,75 GHz frequencies over grass-covered fields. These ob- 
servations were made from an aircraft platform during 
6 flights over a grass watershed in Oklahoma. It was found 
that the coherent scattering component is very important at 
angles near nadir, while the vegetation volume scattering is 
dominant at larger incident angles {> 30 degrees) , The 
results show that least-squares fits to scatterometer data 
can provide reliable estimates of the surface roughness 
parameters, particularly the surface height standard devi- 
ation a. The range of values for a for the 6 flights is 
consistent with a 2 or 3 dB uncertainty in the magnitude of 
the radar response. 



1 . INTRODUCTION 


Recently, there have been many microwave radar measurements 
made of Earth terrain under various surface conditions 
fl-8]. These measurements were usually taken with either 
truck-mounted or airborne scatterometers at various fre- 
quencies, and the resulting data are in the form of radar 
backscattering coefficients a°(6), which are expressed in 
unit of decibels (dB) as a function of incidence angle 6. 
Analysis of the measured radar backscattering coefficients 
can provide valuable information about the surface soil 
moisture content, roughness parameters of the soil surface, 
and the vegetation cover. Theoretical simulations of the 
data can increase our understanding of the manner in which 
microwave radiation is backscattered from multilayer media 
(such as from vegetation-covered soils) . Several theoret- 
ical models [9-14j have been developed to simulate the 
backscattering process, and they provide an excellent means 
of describing the returned radar signals to the scatterom- 
eters. However, the theoretical models are usually compli- 
cated and have many parameters, values of which are 
difficult, if not impossible, to obtain over large natural 
or agricultural fields. From an experimentalist's view- 
point, it would be desirable to have a simple model with a 
parametric description of the scattering media with the 
minimum number of parameters necessary to interpret the 
measured backscattering coefficients. 

The present study was undertaken to develop a simple 
"user's" model for simulating the measured radar back- 
scattering coefficients from vegetation-covered fields in 
conjunction with the data obtained by Jackson et al., 

[1, 2] . The model is based on the theoretical work by Fung 
and Eom [11] , but modified to include the effect of a vege- 
tation canopy. In addition, the Fresnel reflectivity which 
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appears in the model [11] was replaced by calculated soil 
surface reflectivity which was obtained from a radiative 
transfer model [15] , using measured profiles of soil mais- 
turc and temperature. Coherent and incoherent scattering 
are included in the model in addition to the canopy volume 
scattering contributions. There are four parameters; two 
of these parameters specify the condition of the soil sur- 
face (the surface height standard deviation a and the 
correlation length i) , and the other two define the char- 
acteristics of the canopy (the canopy optical thickness, t, 
and the canopy volume scattering factor n) • 

Comparison of the model calculations with the measured radar 
backscatter ing coefficients at the frequencies 1.6 GHz 
(L“band) and 4.75 GHz (C-band) over grass-covered fields 
shows good agreement. The model calculations demonstrate 
that the large magnitudes of the measured a*^(0) at angles 
close to nadir are primarily due to the coherent scattering, 
and that the a°(0) values at large incident angles 
0 > 30*^) can be attributed to vegetation canopy scattering. 
The incoherent scattering contributes to the backscatter ing 
coefficient at all angles for a rough soil surface. 


2. THE MODEL 


I 

( 


A radar pulse reflected from a vegetation-covered soil sur- 
face, is subject to two-way attenuation and scattering by 
the vegetation layer, as shown schematically in Figure 1. 

We assume that the geometrical configuration is symmetric 
with respect to the azimuth angle (j), and that Figure 1 
corresponds to the c{) = 0 case. Backscatter ing occurs at 
®s “ = 0 f where (0, c|)) denote the incident 

direction, and (0 , ^ ) the scattered direction. 

s s 

The backscattering coefficient o°(0) of vegetation-covered 
soils can be written in the form [8, 9], 

O/oi 0//1V . O/Qx -2 t/cos 0 ,,, 

0 (0) = 0 ^ (0) + 0 g (0) e (1) 


where u°(0) is the vegetation backscattering coefficient, 

0 °( 0 ) is the soil backscattering coefficient and x is the 
optical thickness of the vegetation layer. 

Following previous investigations [8, 9], the vegetation 

scattering component cr°{0) c/m be approximated by. 


<(6) = (3. . ,-2t/oos e, 


where ri/ which depends on the canopy water content per 
unit area [8], is a vegetation volume scattering factor. 


For a rough soil surface, the backscattering coefficient 

0 ® ( 0)^consists of two components: the coherent backscatter- 

s 

ing coefficient incoherent backscattering co- 
efficient coherent scattering component 

occurs only in the specular direction (i.e., 0e = 9)^ 

Q 
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Figure 1, Schematic view of the scattering geometry. The 
0 and 0 represent the incident and scattered 

angles, respectively. The thickness of the 

vegetation layer is denoted by d. 
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thus a monostatic radar would not receive any return power 
from the coherent scattering component except for normal 
incidence [14J, However, a radar with finite beamwidth 
antenna pattern can receive both coherent and incoherent 
scatterings, particularly near the nadir direction. 

Thus, one can write the a^(0) in the form, 

S 


a°(6) = 


°ooh(«) 


< 7 ° ( 0 ) 

inc' ' 


(3) 


The coherent scattering of microwave waves from a rough soil 
surface has been investigated by several authors [11, 17]. 

The general form, as a function of both incident and 
scattered angles, can be given by [11, 17], 


°ooh'®' ®s' ♦s* “ ® 

where k is the wave number of the incident wave, <S is the 
Dirac delta function and a is the standard deviation of 
the surface height. The quantities q^, q^, and are 
defined as [11, 14] , 

q = k (sin 0 cos <J) - sin 0 cos <[)) 

X s s 

q = k(sin 0 sin d> - sin 0 sin 4>) .cv 

^y ' s ^s (5) 

q = k (cos 0 + cos 0) 

S' 

a = IR I (cos 0 + cos 0 ) cos((b - d)) , (pp = HH or VV) 
o pp s s 

The delta functions in Equation (4) limit the coherent 
scattering to the specular direction, 0=0 and (j> = <J). 

The magnitude of this coherent scattering along the specular 
direction can be obtained by integrating Equation (4) across 
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the 0„ = e and (jj = (J> direction of the scattered solid angle 
s s 

dQg = dcos0gd(j)g, and the result is approximated by 


^coh<®) " f Ch^^'^'Qs'^^s) ^«s 


= 4ir IR 1^ cos 0 e"^ ® 

PP 


( 6 ) 


2 2 2 
where h = 4k a , and the quantity IR„ I represents the 

Zr 

reflectivity of a smooth surface. 

The coherent backscatter ing coefficient ct^q^( 0) defined in 
Equation (6) corresponds to a monochromatic radar beam. In 
practice, a radar has a transmitting antenna pattern with a 
finite beamwidth, and thus the actual coherent backscatter - 
ing coefficient is distributed according to the transmitting 
antenna pattern G^(0) . Also the received radar power is 
determined by the receiving antenna pattern 0^,(0*) for a 
returned coherent beam located at an angle 0' , measured 
from the center of the antenna beam. 

Assume that the product of the antenna gain patterns is rep- 
resented by the Gaussian form. 


f(0) = G^(0) Gj.(0) = exp 


d(0 - 0 q) 



2i 


(7) 


where 0 is the 3-dB antenna beamwidth and a = 4 Zn2, The 
angle 0^ is the location of the center of the antenna beam. 

Since coherent scattering occurs only in the specular di- 
rection (see Figure 1) , the functional form of receiving 
antenna gain pattern for coherent scattering can be obtained 
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from the one given in Equation (7) by changing (6 - 0^) to 

(6 + 0^), The returned coherent scattering at an angle 

0' =0+0^ can make contributions to the measured back- 
o 

scattering coefficient through an appropriate 'coherent' 
antenna gain pattern, which can be approximated by, 


9c(0) 


Gt(0) 


a(0 - 0 q) 


23' 


exp 


a(0 + 0 q) 


23' 


= exp 


a(0^ + e^)’ 


( 8 ) 


Therefore, the 'measureable ' coherent contribution to the 
backscattering coefficient can be defined by the weighted 
quantity. 


<a 


coh 


( 0 )> 


9c<9) 


coh 


( 0 ) 


(9) 


where a . (0) 
coh' ' 


and g (0) are given by Equations (6) and 
c 


(8) , respectively. 


The incoherent backscattering coefficient a?j^^(0) in Equa- 
tion (3) depends on the statistical properties of a rough 
surface: the surface height standard deviation a and the 

correlation length Z, The latter provides a reference for 
estimating the statistical independence of two points on a 
surface [14], Models for have been developed by 

many authors [11, 14, 16, 17], The one developed by Fung 
and Eom [11] is relatively simple in application and it will 
be further developed in this study to fit the backscattering 
coefficient data. 
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Backscattering coefficient from a rough soil surface depends 
on the surface correlation function of its height distribu- 
tion. Both Gaussian and non-Gaussian correlation functions 
have been used for numerical calculations of backscattering 
coefficient [11/ 14, 23] , For mathematical simplicity, the 
Gaussian form of correlation function has been widely used 
in the computation of backscattering coefficient. In this 
study, we assume that a rough soil surface has a Gaussian 
surface correlation function p(?) and a horizontal correla- 
tion length i, i.e.. 


p(C) = exp(-cV^^) (10) 

where C is a distance between two points on the horizontal 
surface. Then following a similar process as given in ref- 
erence [16] , one can show that the incoherent backscattering 

coefficient PP polarization is given by the 

form [11, 14], 




IRppI^ (1 + sin^e) + Re(RppR*pj^) sin20 


X e 


-h CO 


\ (h cos 
n In 




exp 


n-1 


(kil sin 9) 
n 


( 11 ) 


2 2 2 

where h = 4k a , IR„_I denotes the smooth surface reflec- 

PP 

is the complex conjugate of which is 

a component of the reflectivity. For pp = HH, it can be 
related to R^^^ by the relation [11, 14] , 


tivity, and 


2 sin 0 


V 


2 

COS 0 + W Eg ~ sin 0 


^H1 ' 


( 12 ) 


Ot^Sv ^ Sis 

OF POOSI QUALfTY 

where e„ is the complex dielectric constant of soil. For 
other polarizations., explicit forms of Rppi can be found 
in references [11, 14], 

In applying Equation (11) to the backscatter ing coefficient, 

the quantity a? (0) should be weighted by the antenna gain 
X no 

pattern f(0) given in Equation (7), Similar to <CQoh^®^^ 
in Equation (9), one can define the quantity, 

<a°„o<0)> = £(e) 

as the weighted incoherent backscatter ing coefficient. 

Therefore, the total soil backscatter ing coefficient, 

<cr^(0)> weighted by appropriate antenna gain patterns, can be 
written in the form, 

<a°(9)> = 

I ^ / 2aee \ „ , 

= f(9) exp(- + a?„^(6)l 

where Equations (7) and (8) have been employed in arriving 
at the last step in Equation (14) , 

For comparing with the data, the calculated backscatter ing 
coefficient <a°(0)> from Equation (14) and the vegetation 

5 

backscattering coefficient a°(0) from Equation (2) should be 
averaged over the main beam of the antenna patterns, or more 
precisely, over the illuminated target area bounded by the 
main antenna beam, Allen et al,, [7] have presented the 
detailed description of the geometrical configuration, and 
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our final result of the weighted average total backscatter 
ing coefficient can be written as, 




if 


f(9) 


2aQ0 . 

“coh<«) ®*p(- — 


(15) 


+ cj?nc(0) + <(6) 


tanO do 


where the factor tanO comes from the geometrical configu- 
ration of the illuminated target area [7], and the normali- 
zation factor A is given by; 


A 



f (6) tanO do 


(16) 


Equation (15) will be used to fit the data, as described in 
the next section. 
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3 . THE RESULTS 


The formulas derived in the previous section were used to 
fit the measured backscattering coefficient [1/2] at L and 
C-band frequencies over 4 different grass-covered watersheds 
located near Chickasha, Oklahoma in 1978 and 1980. The data 
for HH polarization, taken with airborne scatterometers from 
an altitude of 300 meter, were given in decibels (dB) at 
incident angles between 5° and 50° at 5° increments. The 
soil texture of the fields is silt loams [18] consisting of 
33% sand, 20% clay, and 45% silt, approximately. Soil mois- 
ture profiles were measured within three depth intervals; 
0-2.5 cm, 2.5-5 cm, and 5-15 cm. Also measured were the 
soil temperatures. These measured soil moistures and tem- 
peratures were employed to calculate the soil dielectric 

. . 2 

constant e and the soil surface reflectivity IR„_I , using 
® PP 

Wilheit's radiative transfer model [15]. These calculated 

reflectivity values were used to fit the data, although 

there would be no significant difference obtained if the 

Fresnel reflectivity was employed. 

Equation (15) was used to fit (by a least-squares criterion) 
the measured backscattering coefficient as a function of 
incident angle 0. There are four adjustable parameters: 
o, if T, and n. The first two (i.e., the surface height 
standard deviation a and the correlation length i) specify 
in a statistical manner the geometrical conditions of the 
soil surface, while the last two (the canopy optical thick- 
ness T and the canopy scattering factor n) describe the 
characteristics of the vegetation canopy, which is assumed 
to form a uniform layer over a soil surface. 

In applying Equation (15) to fit the data, one needs to know 
the 3-dB beamwidth of the antenna pattern. According to 
Wang [19] , the L-band scatterometer had 0 - 9°, and the 
B value for C-band was approximately B - 2.5° [20]. 
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Theoretically, any one of the four parameters can be varied 

to obtain the best fit to the data. However, best fit re- 

sults show the vegetation backscattering coefficient a°(0) 
is relatively small (although it dominates at angles greater 
30^^) , and that one can keep x at a fixed value in fitting 

the data. This is due to the fact that the x value for a 

grass canopy is usually very small and therefore Equa- 
tion (2) essentially reduces to the form a°(0) “ r](l - x/cos 0) . 
Thus good fits to the data can be obtained by varying 
keeping x fixed. 

A previous investigation [21] shows that x is proportional 
to the vegetation canopy water content W (in kg/m ) , and 
that it can be given by the simple form, 

X = cW (17) 

where c is a frequency dependent proportionally constant. 

For L-band, it has been shown that c 0,12 [21] . The x 

values for C-band are 2 to 5 times larger than those for 

L-band. In the present work, the L-band data were fitted by 

keeping x = 0.06, which corresponds to an assumption of 

2 

canopy water content W = 0.5 kg/m , a typical value for 
10-30cm tall grass [22] . For C-band, it was assumed that 
X = 0.12 for all the cases considered. 

With X fixed at the above values, there are only the three 
parameters a, A, and n to be varied to fit the data. Since 
the wave number k always appears in places where a or £ 
occurs in the formulas (see Equations 6 and 11) , it is con- 
venient to take the dimensionless quantities ka and kJl, in- 
stead of a and 1, as the adjustable parameters. 

Comparisons of some typical best fits (at L- and C-bands) to 
the data are shown in Figures 2 to 5 for the four different 
grass fields, in these figures, the solid curves represent 
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Figure 2. Comparison of calculated and measured backscattering coefficients 
The asterisks denote the scatterometer data over site RS , and the 
solid curves represent the calculations obtained with parameters 
listed at the tops. The dashed curve (at the L-band) would resul 
if the coherent scattering was excluded from the calculation. 
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the calculated backscatter ing coefficients (in dB) for the 
HH polarization, and the asterisks denote the data. The 
parameter values used in the calculations are listed at the 
top of each figure, and the soil moisture content (in wt-%) 
of the 0-2,5 cm surface layer is indicated on the lower part 
of each figure, together with the date (month/day/year) of 
the scatterometer measurement. 

Figure 2 displays the results for the site R5, which was 
well managed pasture land [1] , The L-band results are shown 
in part (a) , while those for C-band are given in part (b) , 
The model calculations (the solid curves) in Figure 2 agree 
well with the observations at all incidence angles, except 
at 5°, where the scatterometer data probably contain large 
uncertainties due to instrumental system problems at such 
near nadir angles [20], 

The dashed curve in Figure 2a resulted from excluding the 

coherent backscattering component ct°q^( 0) , This shows that 
the coherent scattering makes large contribution at angles 
near nadir, and its importance can be ignored when the in- 
cidence angle is greater than 15°, 

The ka value used in obtaining the L-band result as shown 
in Figure 2a is smaller than that of C-band (Figure 2b) by a 
factor 3, which is the correct ratio of the k values of 
C-band to L-band if a remains constant, as expected. On 
the other hand, the ratio of the two kil values in Fig- 
ures 2a and 2b does not maintain this 1:3 relationship, thus 
it implies that the correlation length Z, which best de- 
scribes the surface backscatter is still wavelength depend- 
ent, This indicates that, perhaps, additional parameters 
are required to specify the soil surface conditions [7, 11], 

Figure 3 shows the results for site R6 , which was nearly 
identical to R5 in terms of soil condition and vegetation 
cover. The best-fit parameter values listed in Figure 3 are 
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comparable to those in Figure 2 , as expected for two nearly 
identical fields. 

Figures 4 and 5 demonstrate the best-fit results and the 
observations over the two watersheds r 7 and R8, which were 
poorly managed pastures [1,2] , Hydrologically , sites R7 and 
R8 were considered identical. The soils (as shown by the 
soil moistures in Figures 4 and 5) in both sites were very 
dry on the experiment date (September 9, 1980). The calcu- 
lated results as shown by the solid curves in Figures 4 and 
5 agree well with the observations, both in magnitude and 
angular variations. Also, the best-fit ka values (as 
listed at the top of Figures 4 and 5) vary with frequency as 
expected. 

Additional calculations and comparisons with the data are 
given in Appendix A. The parameter values used in the 
calculations of Figures 2 to 4, and those of Appendix A are 
given in Table 1, for both L- and C-bands. The last column 
of Table 1 also lists the soil moisture content (in weight- 
percent, or wt-%) within the 0-2.5 cm surface layer of 
soil. The average value of the parameters are listed in the 
bottom row in Table 1. 

Table 1 shows that the surface parameters for sites R5 and 
R6 have approximately the same numerical values, and that 
those for R7 and R8 are also similar. Generally, the ka 
values for R7 and R8 are larger than those of R5 and r6. 

This is in agreement with the fact that the soil surfaces of 
sites R7 and R8 were more highly eroded and therefore were 
rougher than those of R5 and R6 [1, 2] . 
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Table 1. Best-fit parameters obtained from fits to the 
scatterometer data. Note that the value of 
X = 0.06 i'3 fixed for L-band and x = 0.12 
for C-band. The last column gives the soil 
moisture (SM) within the 0-2.5 cm surface layer. 
The average values of these parameters are 
listed in the bottom row. 


DATE 

SITE 

L-BAND 

C-BAND 

SM 

t 0 - 2 . 5 cm ) 

( WT -%) 

ka 

kS 

*1 

T 

ko 

kS 

t) 

T 

6 / 01 /78 

R 6 

0.09 

2.86 

2.3 X 10 ”^ 

0.06 

- 

- 

- 

- 

21.2 

5 / 12/78 


0.14 

4.62 

3.0 X 10 "^ 

0.06 

0.11 

4.73 

1 . 6 x 10 "^ 

0.12 

15.2 

6 / 30/78 


0.11 

4.40 

3.6 X 10 “^ 

0.06 

0.24 

6.31 

2.8 X 10 “^ 

0.12 

30.1 

6 / 24/80 


0.07 

3.31 

1.4 X 10 “^ 

0.06 

0.21 

4.44 

3.2 X 10 “^ 

0,12 

17.7 

8 / 14/80 


0.08 

3.68 

0.2 X 10 "^ 

0.06 

0.22 

3.74 

0.8 X 10 "^ 

0,12 

2.9 

0 / 09/80 


0.06 

3.24 

0 . 4 x 10 “^ 

0.06 

0.16 

4.20 

1.0 X 10 “^ 

0.12 

7.9 

6 / 01/78 

R 6 

0.12 

4.76 

6.7 X 10 "^ 

0.06 


— 


— 

22.0 

6 / 12/78 


0.22 

5.03 

2.6 X 10 "^ 

0.06 

0.14 

4.89 

1,6 X 10 “^ 

0.12 

16.6 

6 / 30/78 


0.10 

4.16 

6.2 X 10 ”^ 

0.06 

0.26 

4,74 

3.1 X 10 “^ 

0.12 

28.4 

6 / 24/80 


0.10 

4.46 

1 . 7 x 10 “^ 

0.06 

0.26 

4.92 

2.7 X 10 “^ 

0.12 

20.8 

8 / 14/80 


0.13 

4.68 

0.4 X 10 "^ 

0.06 

0.34 

6,14 

0.9 X 10 “^ 

0,12 

1.9 

9 / 09/80 


0.06 

3.90 

0.3 X 10 "^ 

0.06 

0.20 

4.99 

1.0 X 10 “^ 

0.12 

5.6 

6 / 01 /78 

R 7 

0.19 

3.80 

6.1 X 10 “^ 

0.06 

_ 

- 

_ 


19.6 

6 / 12/78 


0.14 

3.36 

4.0 X 10 “^ 

0.06 

0.21 

3.06 

1,0 X 10 “^ 

0.12 

12.7 

6 / 30/78 


0.22 

4,67 

11 , 0 x 10 “^ 

0.06 

0.49 

6.97 

3.8 X 10 “^ 

0.12 

19.9 

6 / 24/80 


0.16 

3.91 

1.6 X 10 “^ 

0.06 

0.40 

6.65 

2.6 X 10 “^ 

0.12 

16.4 

8 / 14/80 


0.16 

4.22 

0.7 X 10 “^ 

0.06 

0.38 

4.71 

0.9 X 10 “^ 

0.12 

1.4 

9 / 09/80 


0.09 

3.76 

0.3 X 10 "^ 

0.06 

0.28 

4.38 

0,7 X 10 “^ 

0,12 

4.4 

6 / 01/78 

R 8 

0.16 

3.20 

14.3 X 10 “^ 

0.06 


_ 


__ 

22.3 

6 / 12/78 


0.14 

3,98 

8.9 X 10 “^ 

0.06 

0.15 

4.20 

2.4 X 10 ~^ 

0.12 

13.0 

6 / 30/78 


0.21 

4,79 

16 . 9 x 10 “^ 

0.06 

0.36 

5.07 

4.4 X 10 “'^ 

0.12 

23.5 

6 / 24/80 


0.42 

6,92 

3.1 X 10 "^ 

0.06 

0.46 

6.05 

3.9 X 10 “^ 

0,12 

16,2 

8 / 14/80 


0.17 

3.70 

1.2 X 10 “^ 

0.06 

0.54 

5.06 

1.1 X 10 “^ 

0.12 

2.2 

9 / 09/80 


0.11 

4.13 

0.6 X 10 “^ 

0.06 

0.40 

5.62 

1.0 X 10 ”^ 

0.12 

7.9 

AVERAGE : 



0.14 

4.15 

4.0 X 10 “^ 

0.06 

0.29 

4.84 

2.1 X 10 “^ 

0.12 
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4. DISCUSSION 


The surface parameter values (for the same site) given in 
Table 1 show some variations from one day to another in the 
same year. It is possible that this is due to the fact the 
surface conditions changed during the period of data acqui- 
sition, However, another possible explanation is that the 
apparent variation in the best fit parameters might be due 
to small errors in absolute calibration of the scatter- 
ometer system from one flight to another. 

To investigate this possibility, we made some studies of the 
backscattering coefficient sensitivity to the parameters for 
L-band. The results are given in Table 2, which lists the 
best-fit parameters values that would result, if the meas- 
ured backscattering coefficients, 0 *^( 6 ), were arbitrarily 
increased by 2 dB and 5 dB, respectively, at all angles. 

Comparison of results in Tables 1 and 2 reveals that in the 
case of a°(0) + 5dB, the ka values (in Table 2) are approxi- 
mately 2 times larger than the corresponding ones in Table 1, and 
that the a°(0) + 2dB case requires about a 50% increment in 
the surface parameter values. Therefore, if the soil sur- 
face conditions of the sites remained the same during the 
data taking period, particularly for the same year, the 
apparent variations in the best-fit parameter value (in 
Table 1) can be attributed to small errors of the scatterom- 
eter system from one flight to tl;e other. Also an assump- 
tion of +2 dB errors in the measured angular distributions 
of the backscattering coefficient would adequately account 
for the range of variation in the parameter values as given 
in Table 1. The n values (C-band) in Table 1 are com- 
parable to results reported by other investigators [ 8 ] at 
the 8.6 GHz frequency for a wheat field, which is assumed to 
be structurally similar to the grass in the pasture. 
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Table 2. Best-fit parameters (L-band) obtained from fits 
the cases of a®(0) + 2 dB, and a®(6)+ 5 dB, 
respectively, as explained in the text. Note 
that T = 0,06 was used in all cases. 


DATE 

SITE 

oO ( 0 l + 2 db 

70 ( 0 ) +6 db 

SM 

{ 0 - 2 . 6 em ) 
(WT %) 

\o 

k£ 

rj 

T 

ko 

kS 

n 

r 

6 / 01/78 

R 6 

0.11 

3.48 

6,7 X lO "^ 

0.06 

0.17 

3.89 

1.3 X 10“2 

0.06 

21.2 

6 / 12/78 


0,24 

6.44 

6.4 X 10 “^ 

0.06 

0.39 

7.63 

1.3 X 10 ^ 

0.06 

15,2 

6 / 30/78 


0.17 

6.26 

6,8 X 10 "^ 

0.06 

0.26 

6.76 

1.4 X 10 "^ 

0.06 

30,1 

6 / 01/78 

R 6 

0.17 

6.73 

12.0 X 10 '"^ 

0,06 

0.27 

6,37 

2.6 X 10 '^^ 

0.06 

22.0 

6 / 12/78 


0.36 

6.38 

4.9 X 10 “^ 

0.06 

0.66 

7,21 

0.9 X 10 “^ 

0.06 

16.6 

6 / 30/78 


0.14 

4.80 

9.3 X 10 “^ 

0.06 

0.22 

6,27 

2.0 X 10 “^ 

0.06 

28.4 

6 / 01/78 

R 7 

0.29 

4.92 

13.9 X 10 “^ 

0.06 

0.41 

5 06 

2.6 X 10 " 2 

0.06 

19.6 

6 / 12/78 


0.18 

3.43 

6.6 X lO "^ 

0,06 

0,26 

3.65 

1,3 X 10 "^ 

0.06 

12.7 

6 / 30/78 


0.36 

6.17 

21.2 X 10 ~^ 

0.06 

0.61 

6.66 

4.1 X 10 ”^ 

0.06 

19.9 

6 / 01/78 

R 8 

0.21 

3.67 

28.1 X 10 ~^ 

0.06 

0.31 

3.06 

6.9 X 10 "^ 

0.06 

22.3 

6 / 12/78 


0.18 

4,22 

14,6 X 10 ~^ 

0.06 

0.26 

4.44 

2.9 X 10 ”^ 

0 06 

13.0 

6 / 30/78 


0.30 

6.69 

29.7 X 10 ”^ 

0.06 

0.44 

6.06 

6.8 X 10 "^ 

0.06 

23.6 
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In addition to the soil surface conditions and the vegeta- 
tion parameters as discussed in the previous section, soil 
moisture content also has a large effect on the backscatter- 
ing coefficient. Figure 6 demonstrates some of the calcu- 
lated results of a®{0) as a function of volumetric soil 
moisture (SM) at the five incideroe angles of 10®, 20®, 30®, 
40®, and 50®, as labeled on the curves. These results 
for L- and C-band were calculated with the average values of 
the best-fit parameters, as given in the bottom row of 
Table 1. 

Figure 6 shows that the backscatter ing coefficient increases 
with SM. Its increment in low 3M region is faster than that 
at high SM (>30%). Below SM » 30%, the a®(0) values in Fig- 
ure 6 vary approximately linearly with SM. Linear regres- 
sion analysis of a®(0) values below SM = 30% was 
performed, and the results are listed in Table 3, which con= 
tains the values of intercept and slope for each incidence 
angle. 

Table 3 shows that the backscatter ing coefficient sensitiv- 
ity to soil moisture, defined as da®(0)/dSM (which equals 
to the slope of the regression line) , decrease as the inci- 
dence angle f increases. At 0 > 40°, Table 3 shows the 
slope values are almost zero (i.e., da°(0)/dSM 0). This 
is due to the fact that the vegetation scattering is the 
primary contribution to the backscatter ing coefficient at 
large angles. The slope values in Table 3 are in good 
agreement with the results reported by Ulaby et al., [6], 

Table 3 also shows that the slopes for both L- and C-band 
are relatively constant up to 20®, but that the intercept 
changes by about 10 dB in L-band and 8 dB in C-band. 

Figure 7 shows the comparison of the calculated backscatter - 
ing coefficients and the scatteroraeter data at 0 = 10®, as a 
function of volumetric soil moisture. The calculated results 
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Table 3 


Linear regression results of backscattering coeffi- 
cient versus volumetric soil moisture (<30%), 
at different incident angles 6. 


0 

(DEGREE) 

L-BAND 

OBAND 

INTERCEPT 

SLOPE 

INTERCEPT 

SLOPE 

B 

-9.6 

0.32 

-6.6 

0.32 

10 

-14.9 

0.31 

-9.5 

0.31 

15 

-17.4 

0.39 

-12.0 

0.28 

20 

-19.9 

0.27 

-14.6 

0.24 

25 

■>22.2 

0.22 

-16.4 

0.16 

30 

-23.6 

0.14 

-17,2 

0.07 

35 

-24.1 

0.06 

-17.4 

0.03 

40 

-24.3 

0.02 

-17.4 

0.01 

45 

-24.3 

0.00 

-17.6 

0.00 

50 

-24.4 

0.00 

-17.6 

0.00 
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Figure 7. Comparison oi calculated and measured backscattering coefficients 
at 0 = lOo. 







{solid curves) are the same as shown in Figure 6, and the 
data (asterisks) are plotted as a function of the soil mois- 
ture within the surface 0-2.5 cm of soil depth. Figure 7 
demonstrates that the agreement between the calculations and 
data is reasonably good within experimental errors. 
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5. CONCLUSIONS 


We have shown that the measured angular distribution of 
backscattering coefficient of vegetation-covered fields can 
be satisfactorily reproduced, using the model developed in 
this study. The model takes into consideration both co- 
herent and incoherent scattering from rough soil surfaces. 

In addition, the vegetation scattering is also included in 
the model and it appears to be dominant at large incident 
angles (i.e., 0 > 30°), The coherent scattering component, 
which is very important at angles near nadir, is introduced 
into the model through the antenna gain pattern with finite 
3-dB beamwidth. The incoherent scattering, which vanishes 
for a smooth soil surface, contributes to the backscattering 
coefficient at all incident angles for rough soil surfaces. 

The ka values obtained from best fits to scatterometer 
data of various sites qualitatively correlate with the de- 
gree of roughness of the soil surfaces. Also, the frequency- 
dependence of the best-fit ka values is in agreement with 
expectation in most cases (Table 1) . This implies that by 
least-squares fit to the scatterometer data, one can obtain 
reliable value of the standard deviation of a rough sur- 
face. However, the kl values do not scale properly with 
wavelength , 

Acknowledgment ; We would like to thank Drs, Bruce Blanchard 
and James Wang for discussions concerning the operations of 
the scatterometers. 
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APPENDIX A 


This appendix presents additional calculations and compar- 
isons of scatterometer data of backscatter ing coefficients 
for HH polarization over grass-covered fields as described 
in the main text of this study. It contains 20 measurements 
at L-band, and 16 measurements at C-band (only L-band meas- 
urements were taken on May 1 , 1978) . 

Each set of calculated and observed results is plotted as a 
function of incidence angle. The parameters used in the 
calculations are listed at the top of each plot. These 
parameters are also listed in in Table 1, The asterisks (*) 
denote the scatterometer data, and the solid curves repre- 
sented the calculated results.. The figures (Figures A-1 
through A-20) are arranged according to the field sites, and 
dates . 
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OKLAHOMA: L-BAND DATA AND BEST-FIT RESULTS 
K51GMA-0.09 KL=2.86 ETA=0.0023 AND TAU=0.0B 



0 20 ^ 60 


INCIDENCE ANGLE (DEGREES) 


Figure A-1 


A- 2 




A- 3 


INCI 














INCIDENCE ANGLE (DEGREES) 















DATA 







OKUHOMA. L-BAhO DATA AW BEST-FIT RESU.TS OKUHOftt. C-BAW DATA AND BEST-FIT 

KSlGm-0.13 KL-A.68 ETA-O.OXM AW TAU-O.OS KSIGMAW.31 KL-5.H ETA -0.0092 am 





INCIDENCE ANGLE (DEGREES) 









BACKSCATTERING COEFFICIENT (DB) 


.iMGtNAL PAGl* SS 
OF POOR QUALITY 


OKLAHOMA: L-BAND DATA AND BEST-FIT RESULTS 
K5 1 GMA-0 . 1 9 KL=3 . 80 ETA=0 . 006 1 AND TAU-0 . 06 



0 20 40 

INCIDENCE ANGLE (DEGREES) 


Q 

ti 

j 


Figure A-11 
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